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Abstract 

The profile of interfaces separating different phases of statistical 
systems is investigated in the framework of renormalized field theory. 
The profile function is calculated analytically in the local potential 
approximation, using the effective potential to two loops. It can be 
interpreted as an intrinsic interfacial profile. The loop corrections 
to the leading tanh-type term turn out to be small. They yield a 
broadening of the interface. 

KEY WORDS: Interfaces, field theory 



1 Introduction 

Many systems of statistical physics can develop interfaces, separating dif- 
ferent phases or substances. These include Bloch walls in ferromagnets, 
binary liquid mixtures, liquid-gas systems and systems of immiscible poly- 
mers [TJ [3 [31 HI [S]. The properties of interfaces are important for the 
behaviour of such systems. They have been investigated both experimen- 
tally and theoretically and are subject of current research, see e.g. [£]. 
For example, if substrates are covered by liquid films in coatings or lubri- 
cants, the stability of the films depends on the interfacial tensions between 
substrate, liquid and vapor. Another example concerns the mechanical 
strength of polymer blends against deformations. The structure of inter- 
faces between domains of the minority component, which are embedded in 
the background of the majority component, is of crucial importance for the 
strength of the blend [7] . 

The profile of an interface characterizes its geometry. It determines 
several other properties such as the interface width and interface tension 
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and therefore plays a central role for the physics of interfaces. Interfacial 
profiles are investigated experimentally through scattering of photons or 
neutrons [5]. Theoretical studies of interfacial profiles include numerical 
Monte Carlo simulations [TUl HU H2] and field theoretic calculations [I] . 

The simplest description of interfaces is based on mean field theory [13] . 
The interfacial profile is described by a continuous order parameter <f>(x) 
representing the difference between the concentrations of the two coexisting 
phases. In the Landau- Ginzburg approach the free energy density is written 
as QU 

£[0] = ^d^ + Vatt). (1) 
In the situation with interfaces the potential is of the double-well type, 

K)W>) = !(0 2 -^ 2 ) 2 - (2) 

The minima of this potential correspond to the two homogeneous phases. 
The mean field correlation length £o is defined through the second moment 
of the correlation function in the mean field approximation. It is given by 
the second derivative of the potential in its minima: 

^ = K(v ))- 1 = -^. (3) 
9o v o 

In mean field theory the interfacial profile is given by minimization of 
the free energy density ([T]) with boundary conditions appropriate for an 
interface. This leads to the differential equation 

A^-V£(0)=O. (4) 

If we choose the interface to be perpendicular to the z-axis, we find the 
typical hyperbolic tangent profile [T5] 

<p{z) = v Q tanh r^ j . (5) 

Its width is proportional to the mean field correlation length £g. The 
parameter zq specifies the location of the interface and will be set to in 
the following. 

Corrections to mean field theory come from fluctuations of the order 
parameter field. They can be calculated systematically in renormalized 
perturbation theory. The fluctuations result in different modifications of 
the mean field result. First of all, higher order corrections change the form 
of the profile from the tanh- function to a different function /. Secondly, 
renormalization of the parameters vq and £o becomes necessary and con- 
sequently the mean field correlation length £ is replaced by the physical 
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correlation length £, which diverges near the critical point with a charac- 
teristic exponent v. A profile of the form 

<Kz)=vf(z/2® (6) 

corresponds to the scaling form proposed by Fisk and Widom [T§]. Fi- 
nally, long-wavelength capillary wave fluctuations lead to the roughening 
phenomenon [17j . It implies a further broadening of the interface such that 
in a three-dimensional system the width of an interface depends logarith- 
mically on the system size and diverges in the limit of an infinite system. 

In this article we consider corrections to mean field theory by means 
of the "local potential approximation" , which is popular in the study of 
non-homogeneous structures in field theory and cosmology. It takes into 
account corrections to the local potential V((j>), whereas contributions in- 
volving higher powers of derivatives are neglected. In this sense, the local 
potential approximation represents the leading contributions in a system- 
atic expansion of the effective action. However, it neglects any non-local 
contribution to the effective action. Therefore it is, in particular, insensi- 
tive to effects related to the system size. The resulting interfacial profile 
is of the scaling form ©. In the local potential approximation the ef- 
fects of capillary waves are not fully taken into account. This is due to 
the fact that in the differential equation determining the interface profile 
only the local value of the profile function but no large-scale properties 
like the system size enter. The profile function therefore describes a kind 
of intrinsic interfacial profile. Including large-scale effects would require 
to build on the full effective potential, which is a non-local functional of 
the profile function. In such a calculation additional contributions to the 
interface profile are to be expected, which depend logarithmically on the 
system size |17j . To define a unique separation of an intrinsic profile from 
its capillary wave contribution in such a situation is a nontrivial subject, 
which has been discussed in various papers, see e.g. (TJ [TSl HJ [TOl EH- I n 
our approach, where broadening due to capillary waves is suppressed, this 
problem does not show up. 

Our calculations are performed in D — 3 physical dimensions in contrast 
to the e-expansions, where D — 4 — e and an extrapolation to e = 1 is nec- 
essary. The three-dimensional approach, in the framework of renormalized 
massive field theory, is based on a systematic expansion in a dimension- 
less coupling [13 [2U] . We calculate the potential to second order in the 
loop expansion. Ultraviolet divergencies are treated by dimensional reg- 
ularization (D — 3 — e), which does not vitiate the fact that the results 
for physical quantities strictly refer to D = 3 dimensions. Renormalization 
of the three-dimensional field theory is performed in the scheme used in 
PT] to two- loop order, employing the results of [HI [23] • The differential 
equation for the interfacial profile is solved analytically in terms of special 
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functions [53]. The resulting profile is compared with results from Monte 
Carlo calculations. 

Previous studies of interfacial profiles in the framework of field theory 
have been made in [23 [23 HE] • In [23 Ml M\ the profile is calculated to 
first order in the e-expansion. This amounts to a one-loop calculation. As 
the expansion is around the four-dimensional case, capillary wave effects are 
not included, and an extrapolation to D — 3 is necessary. The calculation 
of 25J is extended to include the effects of an external field in [35]. In [57] 
the interfacial profile is considered in D = 3 dimensions at one-loop order 
in the presence of an external gravitational field. A functional form of the 
profile is given, including capillary wave effects. The dependence on the 
system size is, however, not considered. 



2 Effective potential 



The effective potential can be calculated in the loop-expansion by standard 
field theoretic methods [23 [33 [3JJ • The starting point is the Hamiltonian 
density 

1 

We write the mean field potential as 



h = -d^d^ + VoW. (?) 



m) -?&r + »? + *!g-»V-tf, (8) 

where the bare mass mo = l/£o is equal to the inverse bare (mean field) 
correlation length. The Hamiltonian is 

H = Jd D xh. (9) 

In the physical situation the number of dimensions is of course D = 3, but 
we shall keep D variable to allow for dimensional regularization. 

The normalized partition function with an external source is defined 
through the functional integral 



Z\3\ = Y Q J-D<f>(x)exp {-H[ct>] + Jd D x 



'xj{x)4>(x) , (10) 



where 



Z = Jv4>(x)exp{~Hm , (11) 

and the usual prefactor (3 — 1/kT is set to 1 by a suitable normalization. 
The free energy is 

W\j]=lnZ\j], (12) 
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and finally the Gibbs potential, which is called effective action in field 
theory, is obtained by a Legendre transformation 



T[4> c ] = W\j}- jd n xj(x)4> c (x), (13) 
where the so-called "classical field" is 

8W\j] 

M*)±-£3- (14) 
6j(x) 

For a constant function <fi c (x) = <f> the effective action yields the effective 
potential V{cj)) via 

-r^= fd D xV((f>). (15) 



More generally, in a derivative expansion the effective action can be repre- 
sented as 

- T[</>] = Jd D x lv(<Kx)) + \z{4>{x)) {d<t>{x)f + ■ ■ ■ | , (16) 

where the dots imply terms with a higher number of derivatives, and 
Z(4>(x)) is a field renormalization factor. 

The effective potential V(4>) is the basic quantity for the discussion of 
spontaneous symmetry breaking and non-homogeneous configurations [32] ■ 
It should be noted that it is not a convex function, as the loop expansion 
also reveals. This is, however, not a failure of the loop expansion, but a 
feature of the effective potential, see e.g. [33 [31] or sec. 2.2 of the textbook 
[35j . The local potential represents the density of the Gibbs potential 
for a small volume element, in which the field is constant, i.e. locally in a 
pure phase. The proof of convexity applies to the extensive thermodynamic 
function, i.e. the total Gibbs potential, and not to the local potential. For a 
proof of convexity in the context of field theory see [3H1 130] • I n a situation, 
where the local potential is not convex, but the Gibbs potential is, the 
system is in a mixed phase and the convexity of the total Gibbs potential 
is the result of the composition of contributions from the coexisting pure 
phases. This is the situation where interfaces are present. Therefore it is the 
local effective potential which is relevant for the local interface structure. 
Another non-homogeneous situation, where the non-convex local potential 
plays a role, is the presence of nucleation bubbles in nucleation theory 
[371 [38]. 

The effective action can be calculated by means of the loop expansion 
in the form 

r[0 = roM+r 1 [fl+r a [0] + ~- ) (17) 
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where the leading term Tq = —TL is equal to the negative Hamiltonian. In 
the corresponding expansion of the effective potential, 



V{4>) = + + + — , (18) 

the leading term Vb is the mean field potential ([5]). 

The first order contribution to the effective action is 

with the differential operator 

^) = -A-^ + f0 2 . (20) 

Setting (f> constant, the first order correction to the effective potential can 
be obtained by Fourier transform as 




(21) 



For small values of the field, \(j>\ < -\/ TO o/.9o, this contribution has an 
imaginary part. Actually, for small values of the field the system is unstable 
due to spinodal decomposition. The boundary value for \<f>\ given above 
represents the mean field spinodal. In the full renormalized theory this will 
be modified or even be replaced by a crossover line [39] . According to [33] 
the imaginary part of the effective potential is to be identified with the 
decay rate of the mixed state. For the calculation of the interfacial profile 
the real part of the effective potential has to be employed. 

The two-loop contribution to the effective potential is obtained through 
the calculation of two Feynman diagrams: 

*<*>— iOO"^""^3"" <2> 



In the two-loop contribution usual ultraviolet divergencies appear. They 
originate from large momentum, short-distance fluctuations, and are not 
related to capillary wave fluctuations. This is reflected in the fact that they 
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are identical to those coming from bulk fluctuations in a situation without 
interfaces. 

The divergencies have to be treated in some regularization and renor- 
malization scheme. We choose to employ dimensional regularization in 
D = 3 — e dimensions. It should be noted that this does not amount to 
an e-expansion, since after renormalization e is sent to zero, whereas in the 
e-expansion one has D — 4 — e and the results have to be extrapolated to 
e = 1. We obtain 

(l)V?M + i + r +fc (i) + .(-i + g)U 



'0 



where 7 = 0.577 ... is Euler's constant. This potential is being used for the 
calculation of the interfacial profile. 



3 The Profile 

In the local potential approximation the effective action is approximated 

by 



r LPA [0] = Jd D x 1^(3)) + J 



2Js (#)r • (24) 

The factor Z3 is the usual field renormalization constant in the broken 
symmetry phase, given by 

Z ^YWY (25) 

with v denoting the minimum of the effective potential. 

The interfacial profile 4>{z) is now obtained in the local potential ap- 
proximation as a solution of the differential equation 

Zz 1 A<t>-V;(4>)=0, (26) 

where V r denotes the real part of the effective potential. For an interface 
perpendicular to the z-axis the appropriate boundary conditions are 

lim 4>{z) = v (27) 

z— too 

lim <j)[z) = -v. (28) 

z — >— 00 

The resulting profile function contains diverging coefficients coming 
from the 1/e term in the second order effective potential (f2"3"|) . These are 
treated by the usual renormalization procedure. We adopt the renormal- 
ization scheme used in |21j to two-loop order. The renormalized mass 
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mii = l/£ is equal to the inverse correlation length £, which in turn is 
defined through the second moment of the correlation function. The field 
4> and its value v in the minimum are renormalized according to 

<I>r{x) = — \= 4>{x) , v R =-^=v. (29) 

V -^3 V -^3 

The renormalized coupling is specified as in |40j through 

9R = ~f- (30) 

In addition we define a dimensionless renormalized coupling according to 

9r / q1 n 

UR = —£-d ■ ( 31 ) 
m R 

Employing the relations given in |22l I23j , the bare quantities mo and go 
are expressed in terms of their renormalized counterparts. This yields an 
expression for the renormalized interfacial profile 4>r{z) depending on the 
parameters tur and ur. In this expression the divergencies are canceled, 
as they should. 

The solution is expanded in the form 

*« w = i/tF + S x,,i) + xAi) + ° ( " y 

^X&- (32) 

Here we have introduced the variable 

z- 2 z . 

The leading order contribution has of course the form of the mean field 
profile 

Xo(z) = tanhz . (34) 
It should, however, be distinguished from the mean field profile, as it con- 
tains the physical correlation length instead of the mean field one. It rep- 
resents a refinement of Landau theory that is consistent with scaling |16j 
and holds near the critical point. 

Due to the fact that the effective potential has a point of non-analyticity, 
the first and second order contributions are also non-analytic. For \z\ 
smaller or larger than artanh(l/-\/3), the functions xi an d X2 take dif- 
ferent forms, which are called xis, X2s and Xih X21, respectively. For the 
first and second order correction we obtain [24] 

12 2 

Yi,(f) = — zsech z -\ — sinhzcoshz tanhz (35) 

v ' 12 9 3 v ; 



(33) 
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11 / + sinh 2 z 

Xil(z) = — £sech 2 z sech 2 z artanh — — 

A v ' 12 12 1 sinhz 

2 2 

H — sinh z cosh z tanh z, 

9 3 



— sinh z \\ — \ + sinh 2 £ ( ^ — ^ sech 2 z I , (36) 



and 



X2s (^) = — ^ z 2 sech 2 z tanh £ 

J 13 17 , 2 . 1 . l2 . 815 l2 , 
+ z — tttt: + tttt: cosh z — — sinh z — sech z 
V 216 216 24 1728 

343 , „ . , „ 2 , 3 . . , . 217 
+ — — cosh z sinh z + — cosh z sinh z — — — tanh z 
648 81 216 

+ i sech 2 z ^\/3 — artanh ar tanh tanh 2j 

1 / 1 3 

H sech 2 z(sinh(4z) — 4z) In ( tanh 2 z 

48 \ \ i ' \2 2 

1 ,2 At. /l-V3tanhz\ T . /l-A/Stanhz 1 

+ — sech z < Li 2 = - Li 2 = 

24 \ ^ l-x/3 J \ 1 + V3 , 

/l + \/3tanhz\ /l + VStanhz' 

"S 1-V3 +Ll2 1 + V3 , 



(37) 



1 / v/ — 5 + sinh 2 z 

— — sech 2 z tanh z artanh 2 , - ; 

144 I smh z 



1 1 sinh 2 



1 , V 2 + smn z \ ( A l2 , 223 , 

H artanh - -{2 — 4 cosh z H sech z 

108 1 sinhz | 1 16 

3. ,2. , . 
+ -zscch ztanhz 



+ (4 cosh z + sech z — 9 sech 3 z) W — ^ + sinh 2 z 
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— h sinh 2 z < z (— 4 cosh z — sech z + 9 sech 3 z) 

2 \l08 v ' 

233 . , . 4 145 , A , . 

smn z cosh z sinh z H seen z tanh z 

648 81 288 

z 2 sech 2 z tanh z 

144 

J 13 17 | 2 . 1 . , 2 . 223 - 

z 1 cosn z sinn z seen z 

V 216 216 24 1728 

^coshzsinhz + lcosh^zsinhz 
648 81 

19 , . 3 , 2 , , . 

— tanh z sech z tanh z 

27 8 

■ sech 2 z ( \/3 — artanh — = I artanh f 



6 V v3/ \\/3tanhz 

1 / 1 3 

H sech 2 z (sinh(4z) - 4z) In h - tanh 2 

48 \ 2 2 

1 2 I / 1 ~ V3 tanh z \ / 1 — \/3 tanh z 

+ — sech z < Li 2 ■= — Li 2 ■= 

24 \ \ 1-VS J \ 1 + V3 / 

/l + \/3tanhz\ /l + y^tanhz' 



(38) 



The contributions of the various orders and the resulting profile function 
are displayed in Fig.[TJ For the renormalized dimensionless coupling we 
have chosen ur — 14.3. In the vicinity of the critical point the coupling 
varies only slowly and is close to the universal fixed point value u* R = 
14.3(1), see [H| for a discussion of numerical and field-theoretical estimates. 

The figure shows that the corrections are small compared to the leading 
term and therefore the interfacial profile at second order can be approxi- 
mately described by a tanh function. It can also be seen that the higher 
order corrections lead to an increase of the interface width. The broad- 
ening of the interface can be specified more quantitatively in terms of a 
broadening factor a by fitting the total profile function to 

X fit (z) = - 1 (^f ) 2 ) tanh(az) , (39) 

where the prefactor is fixed through the asymptotic behaviour of x(z). For 
ur = 14.3 we obtain a value of a = 0.80. 

While the one- loop contribution xi vanishes asymptotically for \z\ — > 
oo, the two- loop contribution approaches a small but finite value, affecting 



10 




Figure 1: The renormalized interfacial profile to second order. The curves 
represent the contributions of the different orders and their sum for the choice 
u R = 14.3. 



the total asymptotic value. It also displays a bump-like behaviour, which 
is due to the non-analyticity at artanh(l/\/3), mentioned above. 

A numerical investigation of the interfacial profile in the three-dimen- 
sional Ising model has been made by Stauffer [9J for a temperature, which 
is one percent below the critical one. The resulting order parameter can 
be described rather well by a tanh function and the numerical fit yields 



= 0.372 tanh(0.1z/a), 



(40) 



where a denotes the lattice spacing. 

In order to compare with our result we have to convert the units by 
using 



1 



1 



m R = - = — 

£ af- 

with /_ = 0.2502(8) 02] and v = 0.630 | 



T -T r 



0.80^ = 0.088-. 
2 a 



(41) 

For T/T c = 0.99 this gives 
(42) 



The broadening factor deviates from the Monte Carlo value by only 
12%. The significance of this coincidence is, however, not clear, due to 



11 



the procedure to identify the interface profile in the Monte Carlo calcula- 
tion. This procedure includes the effects of capillary waves, but neglects 
the effects of bulk fluctuations, whose contributions are central in our cal- 
culation, ft is unclear how the resulting difference in the definitions of the 
interface width affects the comparison. 

4 Conclusion 

We investigated the interfacial profile in the framework of renormalized 
field theory using the local potential approximation. For this purpose the 
effective potential has been calculated in the loop expansion to second 
order. The resulting profile function is obtained analytically. It is of the 
scaling form and can be interpreted as an intrinsic profile. To lowest order 
the profile is of the mean field type. For typical values of the coupling 
constant in the scaling region the higher order corrections are small. They 
imply a broadening of the interface of about 25%. 
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